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ABSTRACT 



We study factorization schemes for parton shower models in hadron-hadron colli- 
sions. As an example, we calculate lepton pair production mediated by a virtual pho- 
ton in quark-anti-quark annihilation, and we compare factorized cross sections ob- 
tained in the conventional MS scheme with those obtained in a factorization scheme 
in which a kinematical constraint due to parton radiation is taken into account. We 
discuss some properties of factorized cross sections. 

§1. Introduction 

Hadron-hadron scattering processes are useful for investigating various prop- 
erties of high-energy reactions, as well as to search for exotic particles. However, 
the leading-logarithmic (LL) order of quantum chromodynamics (QCD) is insuffi- 
cient to evaluate these processes. Thus the next-to-leading logarithmic (NLL) order 
contributions should be taken into account. 

In the evaluation of hadron scattering cross sections, logarithmic contributions 
due to collinear parton production are subtracted from a hard scattering cross sec- 
tion, and they are absorbed into the parton distributions of hadrons. The next-to- 
leading-order (NLO) calculation is necessary in order to remove theoretical ambi- 
guities due to the factorization procedure as well as the choice of the factorization 
scale. 

In actual calculations, Monte Carlo methods are powerful tools for the evalua- 
tion of exclusive processes. The initial state parton radiation can be generated by 
using the parton shower models based on perturbative QCD. One such algorithm is 
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proposed in Ref.l), in which the parton showers are generated according to an algo- 
rithm consisting of a model based on the evolution of momentum distributions with 
the modified minimal subtraction (MS) scheme. In this model, the scaling viola- 
tion of the parton distributions is generated with only information from the splitting 
functions of the parton branching vertices and input distributions at a given low 
energy. It has been found that this method reproduces the scaling violation of the 
flavor singlet parton distributions, up to their normalizations, to an accTiracy of the 
NLL order of QCD. Furthermore, exact NLL-order vertex functions for the decay 
branching processes are taken into account in this model. 

In this paper, we study factorization schemes for the parton shower models in 
hadron-hadron collisions, in which the collinear singularities are subtracted using the 
MS scheme. In order to implement factorization schemes that are appropriate for 
the Monte Carlo methods using the parton shower models to an accuracy of the NLL 
order of QCD, we take into account a kinematical constraint due to virtual partons. 
The momenta of partons in the scattering processes are conserved in this method. 

At NLO accuracy, the processes qq —>■ 'y*g, qg — > 7*g and qg — > 7*5 contribute 
to the hard scattering cross section in actual physical systems. As an example, we 
calculate factorized cross sections for the process qq j*g using our method. Other 
processes can be calculated with similar method presented in this paper. 

In §2, we summarize the calculation of the cross section of the qq — > 7*5 process 
using the dimensional regular ization method. Factorization schemes for the collinear 
singularity are explained in §3. Some properties of the factorized cross sections are 
presented in §4. Section 5 contains a summary and some comments. 



§2. Cross section of the Drell-Yan process 



In this section, we summarize the calculation of the cross section of the Drell-Yan 
lepton-pair production in quark ((j')-anti-quark (q) annihilation. 



Q{Pq) + Q{Pq) 7*(«) + 9{Pg) {P-)+ + gijPg)^ 



(2-1) 



mediated by a photon 7* with virtuality {p- 



q^ = Q'^, where a gluon {g) is 



radiated in the final state. Here, pi {i = q, q, g) and p± denote the momenta of the 
corresponding particles. 

Ignoring the particle masses, the cross section of the Born term. 



Qipq) + q{pq) ^ i*iq) ^ rip-) + i^{p+), 



is given by 



with 



da. 



qq 



NcQ 



3-2er(2-2e) [47r/x2_ 



Q 



2 1 -e 



(2-2) 



(2-3) 
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in 4 — 2e dimensions. Here Nq = 3 and the electric coupling constant of the quark 
is defined by e^afi'^'^, with the mass parameter jj, and the dimensionless coupling a. 
The differential cross section for real gluon radiation is given by 



dfd{-t) 27rs 



'^o(Q',e) 



i-i)i-n) 



with 



I K^^'Jff -t e) 



-P,,(f,e)-2 



where f = Q^/s and 



P,,(f,e) = AW(f) + eP' (f), 



with 



Pi?(.^) = CF^, P;(f) = -CHI - r). 



(2-5) 



(2-6) 



(2-7) 



(2-8) 



Here, Cp = 4/3 and Og are the color factor and the strong coupling constant, re- 
spectively. The Mandelstam variables are defined by 



S = {Pq+ Pqf, i ={Pq-Pgf, U= {Pq - Pgf 



(2-9) 



which satisfy s + i+u = Q'^. 

Integrating over the range < — i < s(l — f) and adding virtual loop contribu- 
tions, we obtain the cross section^) as 



do-qq a,, 



dr 2-iT 



r(l-2e) [47r/x2 



with 



Kq,{r,e) = +4(1 + f^) 



Kqq{f,e), 



2\ /log(l -f] 



+ (_8 + -7r2)5(l-f) 



Here, we have 



(/(r)) ^ = /(f) - (5(1 - f) dyf{y) 
for a function /(f) unregulated at f = 1. 



(2-10) 



(2-11) 



(2-12) 
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§3. Factorization schemes 

In actual Monte-Carlo simulations with the parton shower models, the virtu- 
alities and longitudinal momentum fractions of the partons in the initial state are 
generated according to non-branching probabilities (form factors) and splitting func- 
tions of parton branching processes, respectively. The four-momenta of the partons 
are constructed from these values. 

In order to obtain a finite cross section for the process qq "y*g, we subtract the 
coUinear contributions corresponding to the branching processes q—*qg and q ^ qg 
from the hard scattering cross section. 

The momentum of the quark r^, defined by = Pq—Pg for the branching process 

(l{Pq) Qirq) + g{Pg), (3-1) 

is described by using the momentum fraction Zq, the virtuality r^, and the transverse 
momentum VqT, with pq ■ VqT = Pq ■ Vqx = 0, as 

rq = Pq - Pg = ZqPq + OqPg + VqT, (3-2) 

where a,, = rf^/s for on-shell gluon radiation (p^ = 0). Here, we set Pq=p'q = 0, 
because the relations —Pg, —p^ <^ —r^ are expected in parton shower generation. 
The virtuality of a photon is expressed as 

Q^ = irg+Pgf, (3-3) 
to 0{as) accuracy of QCD. Using Eqs. (3-2) and (3-3), the quantity f is given by 

T = Zq + aq. (3-4) 

According to Eq. (3-4), the subtraction term divided by (To(Q^,e) is defined by 

- ' dzq^r^^Mzq - f - {-rl)/s), (3-5) 



drd{-rl) Jo 'dzqd{-rl) 



with 



_ as 



dzqd{-rl) 27r r(l - e) 



^1 



-rl 



where the splitting function P^\zq,e) depends on the factorization scheme F. 

The conventional MS subtraction scheme^^ corresponds to the approximation 
5{zq — f — {—rg)/s) ~ 6(zq — f) in Eq. (3-5). In this approximation, the collinear 
contribution is subtracted from the hard scattering cross section in the range < 
—t < M^, with f = Zq and —t = — r^. Here, the mass parameter M corresponds 
to a factorization scale for the separation between the initial state radiation and the 
hard scattering process. In the case of the conventional MS scheme, is larger 
than the kinematical boundary for —t, given by 1(1 — f) in 1 — tm < t < 1, where 
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we define tm = M'^ js. Taking the condition = f + {—v^js < 1 into account, the 
colhnear contributions are subtracted within the phase space for the hard scattering 
process, namely < — t < s(l — f). 

We define the integrated contribution as 



i 



dS, 



1 



'i'dTd{-r^) 27rr(l-e) 



(3-7) 



where / denotes the region of phase space for the hard scattering process to be 
considered. Here, sw^^^^^^ is a boundary of the — integration. For example, the 
splitting function with the MS scheme is defined by 



for the branching process presented in Eq. (3-1). 

With zg = f{F = Mi 
gluon radiation is given by 



(3-8) 



With Zq = t{F = MS) and w^^^^^^ = tm, the subtraction term for collinear 



— e 



+ logfM 



(3-9) 



Similarly, with Zg = f + {-r'^)/s {F = MS') and u;['^s'](5) obtain 



log(l - f] 



1 -f 



3 



1 



1 



(l-f)-^7r2(5(l-f) 



(3-10) 



2 (1 - f)_ 

which is valid for 1 — tm < t < 1 {I = S). 

The calculational method employed in the soft gluon region is explained in Ap- 
pendix A. By replacing q with q in the above equations, we obtain the subtraction 
terms for the anti-quark legs, which are the same as those for the quark legs. 



§4. Factorized cross sections 



The factorized cross section in the factorization scheme F is defined by 



df 



?ao(Q^O)i^^l(^)(f), 



with 



'MS), 



f) = Kqg{f,e) - F[fl(^)(e,f) - F^f ^^(e,f). 



for 1 — Tm < t < 1. Here, we obtain 



2 



(4-1) 
(4-2) 



(PW (f))_^ log TM + (-8 + ^TT^^ 6(1 - f)] (4-3) 
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with the MS scheme,^) and 

^PS'](5)^ = Ci. [ - 2(1 + f) log(l - f) 



(1-f) 



+ 2(l-f) 



(4-4) 



with the MS scheme. 

The factorized cross section in the MS scheme depends on the factorization scale 
M. With the MS scheme, the factorized cross section has no M dependence, because 
the subtraction term is integrated over the range 

< -7-2 < s(l 

In order to evaluate the soft gluon contribution, we integrate Eq. (4-1) over 
the range (1 — 77^)5 < < s, with fixed s, where s is generated by the parton 
showers. Here, r]s is a cut-off parameter, satisfying r]s < tm- The energy scale of the 
running coupling constant in this region can be chosen as Q^. In order to simplify 
our analysis, we evaluate the soft gluon contribution with the coupling constant 
Q^sCQ^) — cts(s). The running coupling constant to the accuracy of the NLL order is 
normalized as as (M|) = 0.114 at the Z° boson mass. 

The integrated cross section is given by 

4?(^)(l,l-,.).^#^o(.^o)4f^^(l,l-..), 



27r 



with 



Integrating over the range 1 — % < r < 1, we obtain 



(4-5) 



(4-6) 



r[MS](S) 



(1, 1 - = Cf [4 {5P_fe, 0) - r?,(logr?, - 1) + log^ ,7,} 

-2 j - log(l - 77,) - % + ^ + 2 log % 1 log TM - 8 + ^TT^l (4-7) 



with the MS scheme, and 
.[MS'] (5), 



I^f j^-; (1^ ^_^^^^Cf - 2SP^{vs,0) - {2ns + 3) log vs 



+ log(l-r?,)-8 + -7r^ 



with the MS scheme. Here, the function SP- is defined by 

logf 



SP-{a,b) = / df 



1-f 



We calculate the rjs dependence of the function 



ao{s,0) ZTT 



(4-8) 



(4-9) 



(4-10) 
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by using the expressions appearing in Eqs. (4-7) and (4-8). 

In Fig.l, the ris dependence of the contribution obtained with the MS scheme is 
plotted by the sohd curve. The result obtained with the MS scheme for rjs = tm < 
0.5, which corresponds to the contribution integrated over the range 1 — tm < t < 1, 
is represented by the dash-dotted curve. The rjs dependences found with the two 
schemes become similar, because the log^ ijg term is canceled by the term log rjs log tm 
in the MS scheme, represented by Eq. (4-7). 

The dashed curve plots the rjs dependence obtained with the MS scheme at 
Tm = 0.5. The cross section near the threshold {rjg <C tm) becomes large due to 
the log^?75 term in Eq. (4-7). The error in the calculation is approximated by 
1 < as{{l - r)s)s)/as{s) < 1.05 for < r/^ < tm = 0.5 at = 200 GeV. 



2.5 



MSB(tm= 
MSB(xm= 
MSB' 



0.5) 




0.5 - 



0.01 



0.1 



Fig. 1. The integrated contributions of the cross sections divided by <to(s, 0), which are defined by 
Eq. (4-10), with VS = 200 GeV. The solid curve represents the result obtained with the MS 
scheme. With the MS scheme, the dashed curve and the dash-dotted curve represent the result 
found with tm = 0.5 and that found with tm = rjs, respectively. 

Next, we calculate the factorized cross sections for hard and collinear gluon 
radiation (/ = C) with tm <S 1 — r (see Appendix B). The remnant in the hard 
collinear region integrated over the range < —t < is given by 



^pW(f)log(l-f) + (l-f) 



• (4-11) 



Negative contribiitions remain for f ~ 1 both for the cross section obtained with the 
MS scheme and that obtained with the MS scheme. 

Though the parton showers already include the collinear contributions of the 
g — > qX and g qX processes, as well as those of the q qX and q qX processes, 
at NLL order, ^) we have to add the NLO contributions for the processes qg — > 7*g 
and qg j*q in the hard scattering region. The calculation of the factorized cross 
sections is rather straightforward, because there is no infra-red singularity at f = 1. 
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§5. Summary and comments 



In this paper, wc have studied factorization schemes for hard scattering cross 
sections in hadron-hadron collisions, in which initial state radiation is generated 
by parton shower models. Here, a kinematical constraint due to virtual parton 
contributions is taken into account in the subtraction of collinear singularities. In 
this method, this kinematical constraint guarantees the proper phase space boundary 
for the subtraction terms. 

As an example, we calculated hard scattering cross sections for the process qq — > 
7*5 using the kinematical constraint for the gluon radiation. This method is called 
the MS scheme. We also calculated hard scattering cross section for this process 
without using the kinematical constraint. 

The conventional MS scheme gives a double logarithmic term, which increases 
the cross section in the soft gluon limit, whereas the MS scheme gives only a single 
logarithmic term. The infra-red behavior is rather stable in the case of the factor- 



In the collinear region < — £ < M^, after the collinear singularity is subtracted, 
negative contributions remain for f ~ 1 both for the cross section obtained with the 
MS scheme and that obtained with the MS scheme. Such contributions cannot be 
ignored at NLL-order accuracy. However, event generations employed in the Monte 
Carlo methods according to the negative probability may not be appropriate. We 
will discuss this point in future papers. 

Realizing the matching between the initial parton radiation generated by using 
the parton showers at NLL order and the hard scattering cross sections at NLO 
is not a trivial problem. Further studies of this problem are necessary in order to 
construct consistent algorithms to evaluate higher-order contributions in QCD. 
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Appendix A 



In order to evaluate the quantity 



Jo Jo 



r-{-Tl)/s)CF 



( 




+ 



we calculate 



G[f'K^)(e,f) 
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The function Fjf^^(^) 



(e, f) is given by 



dG, 



[MS'](S) 



99 



df 



Cf 



+ 5(l-f)G[^s'K^)(0) 



.2, /log(l - f) 



-(l-f) + 



+ 

1 1 



7r2(^(l - f) 



2(l-f)+ 3 
which is presented in Eq. (3-10) of the main text. 

Appendix B 



In this appendix, we derive a remnant in the hard colhnear region (I = C) for 
tm ^ 1 — T. The function K^^{T,—i,e) presented in Eq. (2-6) integrated over the 
range < — £ < ( or < —u < M^) is given by 



s Jo 

1 



~ Cf 



ICf 



1 



— e 



logfM(l - f) + 1 - f 



The factorization terms are given by 



F[fl(^)(e,f)= (^ + logfM)Pr(f)^Fr'K^)(e,f 



99 



99 



Thus we obtain the remnant in the hard colhnear region as 



,(C) 



[MS'](C) 



99 



99 

~ Cf 



99 



99 



(r) 



^P,(°)(f)log(l-f) + (l-f) 



which is presented in Eq. (4-11) of the main text. 

The contribution of the hard gluon radiation (I = H) integrated over the range 
M'^ <-i < s{l -t)- M"^ is given by 



= 2Cf 



s{l-f)-M^ 



d{-i)k^^{T-i,e) 



^PS^ir) log ' ! - (1 - f - 2fM)' 

Cf tm 



for f < 1 — 2tm- 
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